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ABSTRACT 

To quantify how rare the bullet-cluster-like high-velocity merging systems 
are in the standard A Cold Dark Matter (CDM) cosmology, we use a large- 
volume (27 h~ 3 Gpc 3 ) cosmological iV-body MICE simulation to calculate the 
distribution of infall velocities of subclusters around massive main clusters. The 
infall- velocity distribution is given at (1 — 3)-R 2 oo of the main cluster (where i? 2 oo 
is similar to the virial radius), and thus it gives the distribution of realistic initial 
velocities of subclusters just before collision. These velocities can be compared 
with the initial velocities used by the non-cosmological hydrodynamical simula- 
tions of 1E0657-56 in the literature. The latest parameter search carried out by 
Mastropietro and Burkert have shown that the initial velocity of 3000 km s _1 at 
about 2i?2oo is required to explain the observed shock velocity, X-ray brightness 
ratio of the main and subcluster, X-ray morphology of the main cluster, and 
displacement of the X-ray peaks from the mass peaks. We show that such a high 
infall velocity at 2i? 2 oo is incompatible with the prediction of a ACDM model: 
the probability of finding 3000 km s _1 in (2 — 3)i?2oo is between 3.3 x 10~ n and 
3.6 x 10~ 9 . A lower velocity, 2000 km s _1 at 2i? 2 oo, is also rare, and moreover, 
Mastropietro and Burkert have shown that such a lower initial velocity does not 
reproduce the X-ray brightness ratio of the main and subcluster or morphology 
of the main cluster. Therefore, we conclude that the existence of 1E0657-56 is 
incompatible with the prediction of a ACDM model, unless a lower infall velocity 
solution for 1E0657-56 with < 1800 km s -1 at 2i? 200 is found. 



Subject headings: cosmology: theory — methods: statistical — large-scale struc- 
ture of universe 
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INTRODUCTION 



The bow shock in the merging cluster 1E0657-57 (also k nown as the "bullet cluster") 
observed by Chandra indicates that the subcluster (found by iBarrena et al.l 120021 ) moving 
through this massive (10 15 /t -1 M M ) main cluster creates a shoc k, and the shock velocity is as 
high as 4700 km s _1 (IMarkevitch et al.ll2002r iMarkevitchl 120061 ) . A significant o ffset between 



the d i stribu tion of X-ray emission and the mass distribution has been observed ( jClowe et al. 



2004 . 120061 ). also indicating a high- velocity merger with gas stripped by ram pressure. 



Several groups have carried out detailed investigations of the physical properties of 



1E0657-57 using non-cosmological hydrodynamical simulatio ns ( )Takizawdl2005l . l2006l ; lMilosavljevic et al 
20071 ; ISpringel fc Farrarll2007t iMastropietro fc Burkertl 120081 ). One of the key input parame- 
ters for all of these simulations is the initial velocity of the subcluster, which is usually given 
at somewhere near the virial radius of the main cluster. 

An interesting question is whe ther the existence o f such a high- velocity merging system 
is expected in a ACDM universe. lHayashi fc White! (120061 ) were the first to calculate th e 
likelihood of subcluster velocities using the Millennium Run simulation (ISpringel et al.ll2005l ). 
As the volume of the Millennium Run simulation is limited to (0.5 h^ 1 Gpc) 3 , there are only 
5 cluster-size halos with Mono > 10 15 hr 1 M n , and 1 cluster with Mono > 2 x 10 15 h^ 1 M m at 



z = 0.28 (close to the redshift of 1E0657-57, z = 0.296). Therefore, lHavashi fc White! (120061 ) 
had to extrapolate their results for M200 > 10 14 h" 1 M & assuming that the likelihood of 
finding the bullet-cluster systems scales with Kub/^ooj where V^ u b is the subcluster velocity 
in the rest frame of the main cluster, and V200 = (GM 2 oo/R2oo) 1 ^ 2 - Here, i? 2 oo is the radius 
within which the mean mass density is 200 times the critical density of the universe, and 
M 200 is the mass enclosed within i? 2 oo- 



While lHayashi fc White! (120061 ) concluded that the existence of 1E0657 -57 is consistent 
with t he standard ACDM cosmology, this conclusion was later challenged by lFarrar fc Rosen 
(120071 ) who showed that, once an updated mass of the main cluster of 1E0657-57 is taken 
into account, the probability of finding 1E0657-57 is as low as 10 -7 . This conclusion still 
relies on the extrapolation of the likelihood derived for M 2 oo > 10 14 h~ x M & . 

As the probability of finding high- velocity merging systems decreases exponentially with 
velocities, an a ccurate determinati o n of t he su bcluster velocity, ra t her t han the shock veloc- 
ity, is crucial. iMilosavljevic et al.l ( 120071 ) and ISpringel &: Farrarl (120071 ) used hydrodynam- 
ical simulations to show that the subc luster velocity can be sig nificantly lower than the 
shock velocity (which is 4700 km s _1 ). iMilosavlievic et al.l (120071) found that the subclus- 
ter velocity can be 4 050 km s -1 , whereas ISpringel fc Farrar ( 2007 ) found that it can be as 



low as 2700 km s 1 . IMastropietro fc Burkertl (120081 ) showed that the subcluster velocity of 
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3100 km s 1 best reproduces the X-ray data of 1E0657-57. 



These varying results are due in part to the varying assumptions about t he initial velocity 



given to the subcluster at the beginning of their hydrodynamical simulations: IMilosavljevic et al. 



( 120071 ) used zero relative velocity between the main cluster and subcluster at the initial sep- 
aration of 4.6 Mpc (which is 2 times -R200 of the main clus ter, 2.3 Mpc). Th e veloc ity is 
about 1600 km s _1 at a separation of 3.5 Mpc (~ l.Sir^ooM Springel fc Farrar ( 2007 ) used 
the initial velocity of 2 57 km s _1 when the separation was 3.37 Mpc (~ 1.5i?2oo)i an d 
Mastropietro fc Burkertl ( 120081 ) explored various initial velocities such as 2057 km s _1 at the 



initial separation of 3.37 Mpc and 2000, 300 , and 5000 km s 1 at the initial separation of 



5 Mpc (~ 2.2i?2oo)- iMastropietro &: Burkertl (120081 ) found that the simulation run with the 



initial velocity of 3000 km s 1 best reproduces the X-ray data. 



In this paper, we demonstrate that the initial velocities used by IMilosavljevic et al. 



but those of 



( 120071 ) and ISpringel fc Farrar! (12007 ) are c onsistent with the prediction of a A CDM model, 



Mastropietro fc Burkert 



(120081 ) at 5 Mpc are not. The simulations of IMilosavljevic et al 



( 120071 ) and ISpringel fc Farrarl (120071) do not reproduce det ails of the X-ray and weak lens- 
ing data of 1E0657-57, and IMastropietro fc Burkertl ( 120081 ) argue that one needs the initial 
velocity of 3000 km s _1 to explain the data. If this is true, the existence of 1E0657-57 is 
incompatible with the prediction of a ACDM model. 



FINDING CLUSTERS OF CLUSTERS IN SIMULATION 



As high- velocity systems are rare, it is crucial to use a large- volume simulation to derive 
the reliable probability distribution. The previous study is somewhat inconclusive due to 
the limited volume of the Millennium Run simulation, (0.5 h^ 1 Gpc) 3 . We calculate the 
probability of finding bullet-like systems using a simulation with substantially larger volume, 
(3 h- 1 Gpc) 3 . 

We use the publicly available simulated dark-matter halo catalogs at z = and 0.5, 
which are constructed from the largest- volume A^-body Marenostrum Institut de Ciencies de 
l'Espai( MICE) simulat ions (jCrocce et al.ll2010l ). They used the publicly available GADGET- 
2 code (jSpringell 120051 ). with the cosmological parameters of fl m = 0.25, £7a = 0.75, Qt = 
0.044, h = 0.7, n s = 0.95, and erg = 0.8. These numbers are consistent with those derived 



from the seven-year data of the Wilkinson Microwave Anisotropy Probe (IKomatsu et al. 



1 Milosavljevic (2010), private communication. All velocities quoted throughout this paper are calculated 
in the rest frame of the main cluster. 
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20101)- 



The MICE simulation that we shall use in this paper has the particle mass of M par = 
23.42 x 10 10 h~ l M and th e linear box size o f L hox = 3072 h' 1 Mpc. The standard friends- 
of-friends (FoF) algorithm ( iDavis et al.lll985l ) with the linking length parameter of b = 0.2 
was employed to find the clust er halos from the di stribution of 2048 3 dark matter particles. 



See iFosalba et al.l ( 120081 ) and ICrocce et al.l ( 120101 ) for a detailed description of the MICE 
simulations and the halo-identification procedure. 

The halos identified in the MICE simulation contain at least 143 iV-body particles. The 
derived halo catalog contains the center-of-mass positions (X) and velocities (V) of halos, 
as well as the number of particles in each halo (N pai ). Note that the number of particles 
in each halo has been co rrected for the known systematic effect of the FoF algorithm, using 



N™ v = iV par (l - iV p -° r - 6 ) (IWarren et al.ll2006l : ICrocce et al.ll2010l ) 



The mass of each halo is calculated as iV n 



The 



.p ar times the mass of each particle, M par . 
mass of halos identified by FoF with the linking length of 0.2 approximately corresponds 
to M 2 oq, i.e., the mass within -R2005 within which the overdensity is 200 times the critical 
density of the universe at a given redshift, M 2 oo = ifpOOpc^Jfi^oo- ft however, known 
that the FoF mass tends to be large r than M 20 o, especially for high- mass clusters which are 
less concentrated (ILukic et al.ll2009l ). As a result, -R200 we quote in this paper may be an 
overestimate. 

The difference betwee n the FoF mass a nd M200 decreases as the number of particles 
per halo, N par , increases (ILukic et al.l 120091 ) . For the main halo masses of our interest, 



A^main > 0.5 x 10 15 h~ x M Q , the average value of N pai is 3355 and 3160 at z = and 0.5, 
respectively. Using this, we estimate that, on average, our i? 2 oo m &y be 10% too large. This 
error is insignificant for our purpose. Moreover, as correcting this error strengthens our 
conclusion by making the probability of finding high-velocity subclusters even smaller, we 
shall ignore the difference between i?2oo and the radius estimated from the FoF mass. 

To find the "clusters of clusters" (i.e., groups of clusters with one massive main cluster 
surrounded by many less massive satellite clusters), we treat each cluster in the catalog as a 
particle and re-apply the FoF algorithm with the linking length of 0.2. This time, the linking 
length of 0.2 means the length of 0.2 times Lb ox /(^ci) 1 ^ 3 , where N c \ is the total number of 
clusters found in the simulation (2.8 and 1.7 million clusters at z = and 0.5, respectively). 
Each cluster of clusters has the "main cluster," or the most massive member of each cluster 
of clusters. All the other clusters are called "satellite clusters" or "subclusters." Table [TJ 
shows the total number of cluster-size halos found in the simulation, the number of clusters 
of clusters having at least two members, and the mean mass of main clusters. For each main 
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cluster, we calculate -R200 from its mass as #200 = [3M 200 /(47r x 200p c (^))] 1 / 3 . Most of the 
satellite clusters are located at r > 2R 2 oo from the main cluster, where r is the distance 
between the main cluster and its satellites. 



3. DERIVING THE INFALL VELOCITY DISTRIBUTION 

Our goal in this paper is to derive the distribution of infall velocities around the main 
clusters. To co mpare with the initial ve l ocities used by the hydrodynamical simulations 



in th e literature (jMilosavljevic et al.l 120071 ; ISpringel fc Farrarl 12007 



2008), we calculat e the infall velocity distribution within (2— 3)fl 2 nn ( Mastropietro &: Burkert 



20081 ). at 1.5^200 flMilosavlievic et al.ll2007l : ISpringel fc Farrarl 120071 ). and at #200- 



Mastropietro fc Burkert 



We define the pairwise velocity of a satellite cluster, V c , as the velocity of the satellite 
relative to that of the main cluster, V c = V main — V sat . When satellite clusters are close 
to the main cluster, V c must be strongly influenced (if not completely determined) by the 
gravitational potential of the main cluster. Thus, V c should depend on the main cluster 

mass, M main . If V c is solely determined by the gravitational potential of the main cluster, 

1/2 

then V c oc M' ai . In reality, however, it is not only the gravity of the main cluster but also 



2005; 


Wane et al. 


2005; 


Wetzel 


2010) 



Figured] shows the distribution of satellite clusters in the log V c — logM main plane (dotted 
line) at z — 0. There is a clear correlation between V c and M main (the larger the M main is, 
the larger the V c becomes), although it is not simply V c oc M main . The dotted line in Figured] 
shows the distribution of all satellite clusters. Next, we shall select the satellite clusters that 
belong to bullet-like systems. We define the bullet-like system as follows: the main cluster 
exerts dominant gravitational force on satellite clusters, and at least one satellite cluster is 
on its way of head-on merging with the main cluster. More specifically, the following three 
criteria are used to select the candidate bullet cluster systems from the clusters of clusters 
at a given z: 



1. Satellite clusters lie between 2.R 2 oo < r < 3i? 2 oo ; an d thus their motion is predominantly 
determined by the gravitational potential of the main cluster, 

2. Satellite clusters are about to undergo nearly head-on collisions with the main cluster: 
|V c -r|/(|V c ||r|) >0.9, and 



3. The mass of satellites is less than or equal to 10% of that of the main cluster, M sat /M n 



< 



1/10, and the main cluster mass is greater than some value, M main > M i 



crit- 
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The third criterion is motivated by the observation of 1E0657-57 indicating that the mass of 
the bullet subcluster is an order-of-magnitude lower tha n that of the massive ma in cluster, 
and the mass of the main cluster is ~ 10 15 h~ l M & (jSpringel fc Farrarl 120071 ). As the 



latest simulation by iMastropietro fc Burkertl ( 120081 ) sh owed that the mass ratio of 6 : 1 best 
reproduces the observed data of 1E0657-56 (also see iNusserl 120081 ) . we have also studied 
the case with M sat /M main < 1/5, finding similar results; thus, our conclusion is insensitive 
to the precise value of the mass ratio. In Figure [2] and [3j we show the distribution of the 
mass ratio, M sat /M main , at z = and 0.5, respectively. As expected, larger- M sat /M main (i.e., 
closer-to-major- merger) collisions are exponentially rare. This makes 1E0657-57 even rarer, 
if the mass ratio is as large as M sat /M main = 1/6. For the rest of the paper, we shall study 
the case of M sat /M ma ; n < 1/10, keeping in mind that 1E0657-57 can be even rarer than our 
study indicates. 

In Figure [lj we show the distribution of satellite clusters satisfying the condition 1 
(dashed line), the conditions 1 and 2 (dot-dashed line), and the conditions 1, 2, and 3 
(solid line). Note that V c of the satellite clusters that satisfy all of the above conditions 

1 /2 

approximately follows V c oc M^ in . This is an expected result, as the satellite clusters in 
this case are basically point masses (nearly) freely falling into the main cluster. We also find 
similar results for z = 0.5. 

In Table [2] and [3], we show the number of bullet-like systems satisfying all of the above 
conditions at z = and 0.5, respectively. At z — 0, about 1 in 3 clusters of clusters with 
M ma x n > 0.7 x 10 15 h^ 1 Mq contains a nearly head-on collision subcluster. At z — 0.5, about 
1 in 5 clusters of clusters with M ma ; n > 0.7 x 10 15 h^ 1 M & contains a nearly head-on collision 
subcluster. Therefore, head-on collision systems are quite common - but, how about their 
infall velocities? 

We calculate the probability density distribution of log V c using the selected bullet-like 
systems (within (2 - 3)i? 2 oo) at z = and 0.5. The results for M main > 0.7 x 10 15 h~ l M Q 
are shown in Figure H] [z — 0) and [5] (z = 0.5). A striking result seen from Figure H] is 
that, of 1135 bullet-like systems shown here for z = 0, none has the infall velocity as high 
as 3000 km s -1 , which is requir ed to explain the X-ray and weak lensing data of 1E0657-56 
(IMastropietro &: Burkertl 120081 ). A lower velocity, 2000 km s _1 , is also rare: none (out of 



1135) within (2 - 3)i? 2 oo has V c > 2000 km s _1 at z = 0. 

We find a similar result for z = 0.5 (Figure [5]): none (out of 78) has the infall velocity 
as high as 3000 km s _1 , and only one has V c > 2000 km s -1 . However, we would need 
better statisti cs (i.e., a bigger simulation) at z = 0.5 to obtain more accurate probability. 



In any case, IMastropietro & Burkertl (120081 ) argued that an infall velocity of 2000 km s 1 is 



not enough to explain the X-ray brightness ratio of the main and subcluster or the X-ray 
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morphology of the main cluster. These results indicate that the existence of 1E0657-56 rules 
out ACDM, unless a lower infall velocity solution for 1E0657-56 is found. 



The significance increases if we lower the minimum main cluster mass. iMastropietro fc Burkert 



(12008h argue that M main ~ 0.5 x 10 15 h' 1 M Q fits the data of 1E0657-56 better. For a lower 
minimum main cluster mass, M ma i n > 0.5 x 10 15 h^ 1 M Q , none out of 2189 bullet-like systems 
at z = has V c > 2000 km s , none out of 186 systems at z — 0.5 has V c > 3000 km s -1 , 
and only one system at z — 0.5 has V c > 2000 km s _1 . 

To examine whether or not the above results depend on the value of the linking length 
parameter, b, of the FoF algorithm used for finding clusters of clusters, we have repeated all 
the analyses by varying the values of b from 0.15 to 0.5. We have found similar results at 
both redshifts, demonstrating that our conclusion is insensitive to the exact values of b used 
for the identification of clusters of clusters with the FoF algorithm. 



To compare with the init ial velocities used by the other simulations ( iMilosavljevic et al. 



20071 ; ISpringel fc Farrarll2007l ). we need to calculate the infall velocity distribution at 1.5i? 2 oo- 
As most of the subclusters are located at r > 2_R 2 oo 5 we have much fewer subclusters in 
(1 — 2)i?2oo- (There are only 191 subclusters within (1 — 2)i?2oo sA z — 0.) To solve this 
problem and keep the good statistics, we shall use the following simple dynamical model to 
convert the results in (2 — 3)i?2oo to those at 1.5i?2oo as well as at -R 2 oo- 

The motion of the subclusters located in (2 — 3)i?2oo is predominantly determined by 
the gravitational potential of the main halo. This is especially true for those in a nearly 
head-on collision course (i.e., nearly a radial orbit); thus, one may treat a selected sub-main 
cluster system as an isolated two-body system. Under this assumption, the pairwise velocity 
at r in < 2i?2oo is given in terms of the velocity at r out > 2i? 2 oo (which is measured from the 
simulation) and the mass of the main halo (which is also measured from the simulation): 

\ t/2/ \ , 2GM main / R 2 Q0 -R200 \ / 1 s 

V- (r in ) = V c (rout) + —5 , (1) 

-"-200 \ r in r out / 

where G = 4.3 x 10 -9 km 2 s~ 2 Mq 1 Mpc is Newton's gravitational constant. 

In Figure H] and [51 we show the probability density distribution of log V c at z = and 
0.5, respectively. The dashed lines show the original distribution for (2 — 3)-R2oo, while the 
dotted and solid lines show the distribution at 1.5i?2oo and i^ no, respectively, computed 



from equation flTJ. We find that the in it ial ve locities used by IMilosavljevic et al.l (120071 ) 



1600 km s 1 ) and ISpringel fc Farrarl (120071 ) (~ 2000 km s x ) are consistent with the 



predictions of a ACDM model: at 1.5i?2oo, 9 (out of 1135) subclusters have V c > 2000 km s _1 
at z — 0, and 16 (out of 117) subclusters have V c > 2000 km s _1 at z = 0.5. However, these 
simulations do not reproduce the details of the X-ray and weak lensing data of 1E0567-56 
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( IMastropietro &; Burkertll2008l ). and thus this agreement does not imply that the existence 



of 1E0567-56 is consistent with ACDM. 

How reliable is this extrapolation of the infall velocity? To check the accuracy of equa- 
tion (pQ), we compare p(V c ) in 2 < r/i? 20 o < 2.4 measured from the simulation and p(V c ) 
at 2.2i?2oo computed from equation JT]). Specifically we use equation (Q]) to calculate the 
velocity at r in = 2.2i? 2 oo from velocities in 2.5i?2oo < r out < 3i?2oo- In Figure EJ, we show the 
measured in 2 < r/_R 200 < 2.4 (dashed line), the predicted p(V c ) at 2.2i? 2 oo (solid line), 
and the original p(V c ) in 2.5 < r/i? 2 oo < 3 (dotted line). We find an excellent agreement 
between the measured and predicted distribution. 



4. DISCUSSION AND CONCLUSION 



Mastropietro fc Burkertl (120081 ) showed that the subcluster initial velocity of 3000 km s _1 
at the separation of 5 Mpc is required to explain the X-ray and weak lensing data of 1E0657- 
56 at z = 0.296. They argued that a lower velocity, 2000 km s _1 , seems excluded because 
it cannot reproduce the observed X-ray brightness ratio of the main and subcluster or the 
X-ray morphology of the main cluster. 

In this paper, we have shown that such a high velocity at 5 Mpc, which is about 2 times 
-R200 of the main cluster, is incompatible with the prediction of a ACDM model. Using the 
results at z = and M main > 0.7 x 10 15 hr l M , ACDM is excluded by more than 99.91% 
confidence level (none out of 1135 subclusters has V c > 2000 km s _1 in 2 < r/_R 2 oo < 3). 
For a lower minimum main cluster mass, M ma ; n > 0.5 x 10 15 h" 1 M , ACDM is excluded by 
more than 99.95% confidence level (none out of 2189 subclusters has V c > 2000 km s _1 in 
2 < r/R 2Q0 < 3). 

The results at z — 0.5 are not yet fully conclusive due to the limited statistics: none 
out of 78 subclusters has V c > 3000 km s _1 in 2 < r/i? 200 < 3, while there is one subcluster 
with V c > 2000 km s" 1 in 2 < r/R 200 < 3. For M main > 0.5 x 10 15 h' 1 M , none out of 186 
subclusters has V c > 3000 km s _1 , while there is one subcluster with V c > 2000 km s _1 . 

While these confidence levels are directly measured from the simulation, one can es- 
timate the probability better by fitting the probability density, p(\ogV c ), to a Gaussian 
distribution as 



p(log V e ) = ^=L=exp 

\/27T0-~ 



(logK-z/) 2 
2a, 2 



(2) 



where V c is in units of km s 1 and v and a v are the two fitting parameters. The best-fit 
values of the two parameters for z = and 0.5 are (v, a v ) = (3.02,0.07) and (3.13,0.06), 
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respectively. The mean velocity at z = is smaller than that at z — 0.5 by a factor of 
2Q3.i3-3.02 _ ^ 29. This may be understood as the effect of A slowing down the structure 
formation at z < 0.5. 

Generally, one has to be careful about this approach, as we are probing the tail of the 
distribution, where the above fits may not be accurate. Using the above Gaussian fits, we 



find P(> 3000 km s~ 
find P(> 2000 km s" 



3.3 x lO" 11 and 3.6 x 10" 9 at z = 
2.9 x 10~ 5 and 2.2 x 10~ 3 at z 



and 0.5, respectively. We also 
= and 0.5, respectively. These 



numbers pose a serious challenge to ACDM, unless one finds a lower velocity solution for 
1E0657-56. Here, a "lower velocity" may be somewhere between V c < 1500 and 1800 km s" 1 
at r ~ 2i?2oo ; which give 1% probabilities at z = and z = 0.5, respectively. 

The bullet cluster 1E0657-56 is not the only site of violent clu ster mergers. For ex ample, 
there are A520 dMarkevitch et al.ll2005h and MACS J0025.4-1222 faradac et aHl2008h . Also, 
high-resolution mapping observations of the Sunyaev- Zerdovich (SZ) effect have revealed a 



violent merger event in R X J1347-1145 at z = 0.45 ( iKomatsu et al.l 1200 



2004: 



Mason et al.l 120091 ). which are confirmed by X-ray observations (lAllen et al. 



Kit avama et al. 



2002 



Ota et al.l 120081 ) . The shock yelocity inferred from the SZ effect and the X-ray data of 



RX J1347-1145 is 4600 km s -1 (Kitavama et al.ll2004f) . which is similar to the shock velocity 
observed in 1E0657-56 ( IMarkevitchl 120061 ) . The lack of structure in the redshift distribution 
of member galaxies of RX J134 7-1145 suggest s that the geometry of the merger of this 
cluster is also closer to edge-on ( ILu et al.l 120101 ). However, the lack of a bow shock in the 
Chandra image may suggest that it is not quite as edge-on as 1E0657-56. In any case, it 
seems plausible that there may be more clusters like 1E0657-56 in our universe. This too 
may present a challenge to ACDM. 

Since the volume of the MICE simulation is close to the Hubble volumejfl our results can 
be compared directly with observations, provided that detailed follow-up observations are 
available for us to calculate the shock velocity, gas distribution, and dark matter distribution. 
These three observations would then enable us to estimate the mass ratio and initial velocity 
of the collision which, in turn, can be compared to the probability distribution we have 
derived in this paper. Note also that the probabilities obtained in our work are the conditional 
ones. That is, the probability for which a fitting formula is provided is the probability of the 
velocity of bullet-systems that are nearly head-on, with 1:10 or more mass ratio, and with 
15 /i -1 M . If we computed the probability of finding high-velocity bullet 



M main > 0.7 x 10 



2 For example, the comoving volume available from z = to z — 1 over the full sky is 54 h~ 3 Gpc 3 , which 
is only twice as large as the volume of the MICE simulation. The comoving volume out to z = 3 is still 
396 h~ 3 Gpc 3 , which is nowhere near enough to overcome the probability of 10~ 9 . 
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systems among all clusters from the simulation, then the probability would be even smaller 
than those estimated above. Such a conditional probability is relevant to the observation, if 
we have sufficient amount of data for estimating the mass ratio and the initial velocity, as 
mentioned above. Note that we have precisely such data for 1E0657-56. 

An interesting question that we have not addressed in this paper is how many high- 
velocity bullet systems are expected for flux-limited galaxy cluster surveys, such as the South 
Pole Telescope and eROSITA (extended ROentgen Survey with an Imaging Telescope Array). 
To calculate, e.g., dN^ n n et /dz, one needs the light-cone output of the MICE simulation. While 
we have not investigated this, we expect two major light-cone effects on the infall velocity 
distribution. First, the infall velocities at high z's should be larger since the effect of A has 
yet to kick in at high z's, which we have already demonstrated here by comparing the mean 
infall velocity at z = and 0.5. Second, the massive bullet systems with mass greater than 
10 15 /i~ 1 M are very rare at higher z's. The first effect will make the high-velocity system 
more common, while the second effect will make the high-velocity system less common. In 
order to quantify the net effect, one needs the light-cone output. However, the light-cone 
effect alone would not be able to reconcile the existence of 1E0657-56 with the prediction of 
ACDM. 
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13 14 15 

l°g[M maln /(h- 1 M )] 

Fig. 1. — Distribution of satellite clusters in the K-M main plane. The dotted line shows all the 
satellite clusters found in the simulation at z = 0; the dashed line shows those lying between 
2-R200 < r < 3i?2oo from the main cluster; the dot-dashed line shows those lying between 
2-R200 < r < 3i?2oo and about to undergo nearly head-on collisions with | cos#| > 0.9; and 
the solid line shows those lying between 2i? 2 oo < r < 3i?200) about to undergo nearly head-on 
collisions, and having small masses compared to the main cluster mass, M sat /M main < 1/10, 
where M main > 0.7 x 10 15 h' 1 M . 
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Fig. 2. — Distribution of the sub-main cluster mass ratio, M sat /M main , at z — 0. The solid 
and dashed lines show the distribution for the main cluster masses of M main > 0.7 and 
0.5 x 10 15 hr x Mq, respectively. The distribution is normalized to unity when integrated, 
Jq p(x)dx = 1. 
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Fig. 3. — Same as Figured but for z = 0.5. 
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Fig. 4. — Probability density distribution of the infall velocities, log V c , of the bullet-cluster- 
like systems at z = 0. The main cluster masses are M main > 0.7 x 10 15 h~ l M & , for which 
there are 1135 bullet-like systems in the simulation at z — 0. The dashed line shows the 
distribution of log V c within 2 < r/i?2oo < 3 measured from the simulation. This distribution 



shows that the initial velocities used by lMastropietro fc Burkertl (120081 ). V c > 2000 km s 1 at 



2.2i?200) are incompatible with the prediction of a ACDM model: none (out of 1135 eligible 
samples) has the velocity as high as V c > 2000 km s -1 in 2 < r/i? 2 oo < 3. The dotted and 
solid lines show the distribution of V c at 1.5i?2oo and -R200, respectively, which are obtained 
by converting the dashed line using equation ((Tj). We also show a Gaussian fit to the dashed 
line, which is given by equation fl2}. Note that 10 3 2 = 1585, 10 3 - 3 = 1995, and 10 3 - 4 = 2512. 
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Fig. 5. — Same as Figure HJ but for z = 0.5. The main cluster masses are M main > 0.7 x 
10 15 h^ 1 Mq, for which there are 177 bullet-like systems in the simulation at z = 0.5. 
None has the velocity as high as ^ = 3000 km s _1 in 2 < r/i?2oo — 3, while there is one 
subcluster with the velocity of V c > 2000 km s _1 in 2 < r/i? 2 oo < 3. (This subcluster 
has V c = 2049 km s -1 .) We also show a Gaussian fit to the dashed line, which is given by 
equation (J2J. Note that 10 3 2 = 1585, 10 3 3 = 1995, and 10 3 - 4 = 2512. 
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Fig. 6. — Testing equation (p]). The dashed line shows the distribution of logVc in 2 < 
r /-R2oo < 2.4 measured from the simulation, while the solid line shows the distribution of 
log V c at r in = 2.2i? 20 o calculated from the measured distribution in 2.5 < r out /i? 2 oo < 3 
(dotted line) and equation (Tj[|). 
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Table 1. 



z 


the number of 


the number of 


the mean mass of 




clusters 


clusters of clusters a 


main clusters 





2.8 million 


0.29 million 


1.3 x 10 14 h- 1 M e 


0.5 


1.7 million 


0.20 million 


1.1 x 10 14 h- 1 M Q 



a A "cluster of clusters" is a group of cluster-size halos identified by the 
FoF algorithm. A useful picture is a massive cluster surrounded by many 
less massive clusters. 
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Table 2. 



-^main 


the number of 


the number of 


the number of 


[io 15 ^-^©] 


clusters of clusters 


bullet-like systems 51 


bullet-like systems' 3 




at z = 


at z = 


at z = 


> 0.5 


8523 


2189 


3093 


> 0.7 


3135 


1135 


1402 


> 1 


911 


351 


391 



a For M sat /M main < 1/10. A "bullet-like system" is defined as a nearly head-on collision 
system satisfying all of the conditions (1, 2, and 3) given in Section 3. 

Tor M sat /M main < 1/5. 
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Table 3. 



-^"main 


the number of 


the number of 


the number of 




clusters of clusters 


bullet-like systems 3. 


bullet-like systems b 




at z — 0.5 


at z — 0.5 


at z — 0.5 


> 0.5 


3108 


186 


240 


> 0.7 


800 


78 


93 


> 1 


138 


27 


32 



a For M sat /M main < 1/10. 
b For M sat /M main < 1/5. 



